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Abstract 

Diffuse scattering is usually associated with some disorder in the analyzed 
material. Different kinds of disorder may produce different diffuse scattering 
- or not. In this letter, we demonstrate some aspects of the variety of diffuse 
scattering that occurs even in very simple examples, and how unawareness may 
lead astray. 



1 Introduction 

For a long time, the diffuse part in diffraction spectra played only a minor role in 
crystallography. This was mainly due to experimental restrictions but also to a lack 
of theoretical studies. In this letter, we investigate some simple models with disorder, 
both deterministic and random, where the diffuse background of the diffraction must 
not be neglected in the structural analysis. 

We first start with an introduction to the language of mathematical diffraction 
theory, which is necessary to get a full understanding of the spectrum. The following 
two examples have exactly the same diffraction pattern albeit their disorder is of com- 
pletely different type. We then discuss models with singular continuous diffraction, 
a type of spectrum that is not present in classical crystallography but may appear 
in random structures with long range order. This is followed by the Ising lattice gas 
where the Bragg spectrum alone does not reflect the correct symmetry. Last, we 
compare square ice with an alternative model that can only be distinguished by the 
diffuse scattering. 
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2 Diffraction theory 



Let us recall some notions of diffraction theory || ||, adapted to our purposes using 
the measure theoretic approach || . On the assumption of kinematical diffraction in 
the Fraunhofer picture, the diffracted intensity per atom %, is given by the Fourier 
transform of the autocorrelation 7^. Let the unit point measure (Dirac distribution) 
8 X idealize an atom at position x and A denote the point set of all possible atomic 
positions. Then the atomic structure is described by the so-called weighted Dirac 
comb 

= y^ j w(x)6 x , (1) 

where w(x) e {0, 1} indicates whether the point x is actually occupied or not. More 
generally, w(x) can be a bounded complex weight function. The autocorrelation 
(Patterson function) 7^, assuming it is well defined, is given by 

7* = (2) 
zeA 

with A = A — A being the set of difference vectors which we assume discrete and 
closed. The autocorrelation coefficient v{z) can be calculated according to 

v(z) = lim 1 V w(y)w(z-y), (3) 
v ; yeA R 

z—yeA 

where B R is the ball of radius R around the origin, A R — A n B R and — complex 
conjugation. If w(x) G {0,1}, v{z) is the frequency of having two scatterers at 
distance z. 

The Fourier transform of a function is 

J(jfe) = f e - 2mk - x <j)(x) dx, (4) 

and we use the standard theory of tempered distributions from here, see || for our 
conventions. 

The diffraction image, is a positive measure that tells us how much intensity 
is scattered into a given volume element. So, by the general decomposition theorem 



T6| , Ch. 1.4], it admits a unique decomposition with respect to Lebesgue's measure 
into three parts, % = (%) pp + (%) sc + (%)ac- The pure point part (%) pp consists 
of the Bragg peaks, the absolutely continuous (7 w ) ac corresponds to the usual diffuse 
background that can often be described as a continuous function, otherwise as an 
/^-function. The singular continuous component (7u;)sc lies somewhere in between, 
i.e. it is neither a smooth function nor "as singular" as a Bragg peak, for a precise 



definition see |T6| 
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3 Random versus deterministic disorder 



In the following, we present two one-dimensional examples with completely different 
type of (dis)order but displaying the same diffraction spectrum, i.e. they are homo- 
metric. 

The first model is a Bernoulli system on Z. The structure is given by the stochastic 
Dirac comb u B = Ylmez V B { m )^m where r] B (m) is a family of random variables that 
take the values hi and h 2 with probabilities p\ and p 2 . For such a lattice gas, one can 
prove || H that its diffraction measure almost surely consists of a uniform pure point 
and a constant absolutely continuous part, 

with mean (h) = p±hi + p 2 h 2 and variance Vax(h) = (\h\ 2 ) — \(h)\ 2 . For hi = 1, 
h 2 = and p l = p 2 = 1/2, we have \(h}\ 2 = Var(h) = 1/4. 

Let us compare this with the Rudin-Shapiro sequence, where the atoms are dis- 
tributed deterministically. We construct u RS = J2 m &z V { m )$m according to the 
rule 

rj RS (0) = 1, V RS (-1) = 0, 
v RS (4m) = r] RS (Am + 1) = r] RS {m), (6) 

V RS (4m + I) = i (l - (_i)-+<+^ s M ) , Z = 2, 3 

Alternatively, this sequence may be defined via a substitution rule on a four letter 
alphabet [H]. For the initiate, we add that we use the square of the traditional substi- 
tution to get the fixed point of the bi-infinite sequence. The diffraction spectrum can 
be calculated rigorously |D| pp. 165-168]. With our choice of the scattering strenghts, 
this results in 

mgZ 

Amazingly, though the structure is completely deterministic, its two-point correlations 
are destroyed systematically so that only a constant diffuse background remains in 

Therefore, with hi — 1 and h 2 = 0, we have 7<f = ^y RS , a distinction of the 
original structure on the basis of the diffraction spectrum is impossible. Even for 
finite systems, the slight difference in the fluctuations of the background is almost 
undetectable. Since the Fourier transform is unique, this means that the Rudin- 
Shapiro sequence is homometric to the above Bernoulli chain (a statement which is 
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true with probability one) . The remarkable feature of this example is the inclusion of 
a diffuse background, see |17| for other cases of homometry. 

Note that this example explores the full entropy range: the Bernoulli case has 
entropy log 2, the maximal value for a binary system, while Rudin-Shapiro has en- 
tropy 0. It is clear that one can find other examples with the same diffraction image 
but entropy between these extremal values. So, a resolution of the corresponding 
reconstruction problem, unless extra information is available, needs an optimization 
approach, e.g. by choosing the structure which maximizes the configuration entropy. 



4 Singular continuous spectra 

There is no simple characterization of the singular continuous part of a diffraction 
spectrum, since this term covers the whole range (in the sense of tempered distribu- 
tions) between a continuous function, resp. /^-function, and a Dirac distribution. As 
indication one may use the scaling behaviour || of the intensities with the system 
size N: 7 PP ~ N, ^y ac ~ const, and ■y sc ~ N a , < a < 1, but this recipe is not 



correct in general |10|. A classical example where the scaling argument yields the 
right answer is the Thue-Morse sequence |L5], |[. Singular continuous spectra usually 
appear in structures where the constructive interference that is responsible for the 
Bragg peaks is impaired by some randomness. Nonetheless, some long range order is 
still strong enough to prevent the peaks from becoming completely diffuse. Note that 
singular "peaks" can never be isolated - a property that further complicates their 
analysis. 

As an illustration in one dimension, let us consider a "tiling" of the line without 
gaps and overlaps which consists of two prototiles (intervals) of length 1 and r = 
(1 + \^E)/2. The unit scatterers are located on the left endpoints of the intervals. By 
placing the intervals randomly with arbitrary (but fixed) probabilities, the resulting 
spectrum is absolutely continuous except for the trivial Bragg peak at k = 0, see || 
Thm. 2] . On the other hand, arranging these intervals as a Fibonacci sequence leads 
to the well-known quasicrystal with purely discrete diffraction spectrum, cf ||. 

With the same tiles, one can now also construct a "structure intermediate between 
quasiperiodic and random". This is achieved using the so-called circle sequence ||. 
The positions x n of the atoms are given by 

x n - x n _i = 1 + £l[ 0!/3) (raa) (8) 

where < j3 < 1, l[o,/3) is the characteristic function of the interval [0,(3) and xo, 
£ G K. are parameters; assume xo = for simplicity. In order to have the same tile 
lengths and density d = 1/(1 + /?£) as in the Fibonacci case, we choose £ = r — 1 and 
(3 = 2 — r. As was shown by Hof ||, the circle sequence has purely singular continuous 



4 



diffraction spectrum apart from the trivial Bragg peak for every irrational £, every (3 
and generic ("most") a. Nevertheless, the assessment in each explicit case is open. 

In two dimensions, a straightforward generalization of the above mentioned Thue- 
Morse chain is given by the next example. It is constructed by using a two-dimensional 
substitution rule |[| on a two-letter alphabet: 



b a 
a b 



p : a i — , : b 
or, using two elementary cells with different atoms 



a b 
b a' 



(9) 



□ 



o • 
• o 



Figure 1: Construction of the 2D Thue- Morse pattern. 

The resulting structure is a simple Cartesian product of Thue-Morse chains and 
can be treated using the method indicated in || Sec. 4]. The diffracted intensity is 



a product measure 7^ u = 7^ • 7™ of the corresponding one-dimensional intensities 

7, 



■TM 



(compare ||). We assign scattering strenghts of +1 (full circles) and —1 (open 
circles) to the atoms. Taking values 1 and would just lead to an additional, trivial 
Bragg part on Z 2 . The remaining diffracted intensity is a product of singular contin- 

and thus itself purely singular continuous (Fig. 



uous measures 
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We have to 

emphasize that this singular continuous behaviour is by no means recognizable from 
the shape of the peaks and the decrease of their wings. The most likely fit is given 
by Gaussians, but this is also the case for true Bragg peaks. 

The explicit scaling of the peaks with the system size can be obtained by applying 
the formulae of || Sec. 4.3] to our product structure. 



5 The Ising lattice gas and its symmetry 

It's worth taking a glance at the well-known 2D Ising model where, in contrast to the 
previous examples, we have some interaction between the atoms. In the lattice gas 
formulation, we place atoms and holes (scattering strengths 1 and 0) on the square 
lattice subjected to coupling constants Ki = Ji/k B T > in x- and y-direction. There 
is a phase transition at (sinh(2i^i) sinh^i^)) -1 = 1- One can show || that the pure 
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Figure 2: Diffraction image of the 2D Thue-Morse product structure (top). Cut along 
the line k — g (bottom). The system size for the numerical calculations was 512 x 512. 



point part of the diffraction spectrum is fourfold symmetric for all Ki and has the 
form 

(%) PP = l% k&25k ; T T -J; do) 

where p is the ensemble average of the number of scatterers per unit volume, 1/2 < 
P<1- 

The correct (two-fold) symmetry is visible only from the diffuse part of the spec- 
trum, the Bragg part alone is misleading. Using the asymptotic behaviour of the 
correlation functions, the background can be shown to be an absolutely continuous 
measure B, which can be represented by a continuous function. The background 
concentrates around the Bragg peaks, as expected for an attractive interaction. 
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Figure 3: Ising lattice gas above T c with K\ = 0.35 and K 2 = 0.1. 

6 The ice model 

Square ice is another example for which the relevant structural information of the 
diffraction image is only contained in the diffuse background. In this model, oxygen 
atoms are sitting on the vertices of a square lattice, see Fig. |[ The hydrogens are 
placed on the edges with a distance of 1/3 to the next O according to the following 



ice (or Bernal- Fowler) rules ([14] and references therein) 



1. There are two hydrogens adjacent to each oxygen. 

2. There is only one hydrogen per bond. 



o -o o -o 

•o 6- o 6 

6- 6 -o -o 
• • 

•6 *o *o 6 

Figure 4: Typical square ice configuration 



By using the equivalence to a special case of the six-vertex model, Lieb |12| solved 
the ice model and calculated its residual entropy to be |log|. Although calorimetric 



measurements on three-dimensional ice confirmed the predicted value of the entropy 
very well, it was almost impossible to detect the hydrogen disorder in X-ray scatter- 



ing experiments. Neutron studies are more promising, compare [JTJ where a highly 
disordered hydrogen structure was shown in deuterated cubic ice Ic. 

Apart from the weak scattering power of the hydrogens, another more fundamental 
problem arises. Consider the case of a Bernoulli system where we ignore the first ice 
rule, so that the hydrogens are distributed with equal probability independently of 
the others. 

By assigning scattering strengths ho and /ih to the atoms, a straightforward cal- 
culation yields the autocorrelation of the Bernoulli model: 

7- = [ h o + h l + -f 22 <W 



x=±\ 



x=±\ 



Thus, the Bragg part is given by {k = (hi, /C2)) 

\ ( ( 2irki 27r/c 2 \\ 2 e 

M P p = 2^ [ h ° + hn \ cos + cos— — J 1 5 k , (12) 

and the absolutely continuous background can be represented by the continuous func- 
tion 

(U\ 1,2/. 2^1 , • 2 * k 2\ 

a( k ) = M sin — + sm — I. (13) 

But due to the same averaged occupation probabilities of the hydrogens in the ice 
and the Bernoulli model, their pure point parts are indistinguishable. In other words, 
Eq. fll2] ) is valid for the case of square ice, too. The background is different, of course, 
but very weak, in particular because, in reality, /io > /ih. 

Fig. [5] displays the diffraction image of square ice. To focus on the contribution of 
the hydrogens, we set ho = 0. Positive ho would result in an increase of intensity of the 
Bragg peaks, except for the case of ho = hn where we have extinction at (k\, fc 2 ) = 



(1,1), (1,2), (2,1) and (2,2), according to Eq. fll2"|) . The continuous background, 
however, depends only on /i H , as in the Bernoulli case. 

We have no explicit formula for the background in the ice model, since the correla- 
tion functions are not known, but we can calculate the diffraction image numerically 
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Figure 5: Elementary cell of the diffraction pattern for square ice (ho = and hn = 1). 
Only the weak (even for vanishing ho\) background differs from the Bernoulli system, 
compare Eq. (O). 



(FFT) from a pattern of a Monte Carlo simulation, see |13|, Ch. 7] for the method. 
The decay of the correlations suggests that the background is an absolutely continuous 
measure and there is no singular continuous component present. Note that the inter- 
pretation of the diffuse part is more involved here - it needs to include multi-particle 
effective interactions. 
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